Finite-size and confinement effects in spin-polarized trapped Fermi gases 
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We calculate the energy of a single fermion interacting resonantly with a Fermi sea of different- 
species fermions in anisotropic traps, and show that finite particle numbers and the trap geometry 
impact the phase structure and the critical polarization. Our findings contribute to understanding 
some experimental discrepancies in spin-polarized Fermi gases as finite-size and confinement effects. 
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Experiments with spin-polarized Fermi gases [J 0, 0, 0, 
HI, S 0] enable a unique exploration of superfluidity and 
universal properties in strongly-interacting asymmetric 
Fermi systems. There are very exciting experimental re- 
sults of the MIT [H,[i,H and Rice University groups, 
however with differences in the observed phase structure 
and the critical polariztion. In this Letter, we provide 
a first microscopic explanation of the MIT-Rice differ- 
ences: The particle number and the trap geometry affect 
the interaction energy in the normal polarized phase and 
this impacts the limit of superfluidity in traps. 

The MIT experiment [3, 0, 0] observed phase separa- 
tion in the trap, with equal densities in the core, sur- 
rounded by a partially-polarized shell and an outer re- 
gion of normal majority fermions. The study of vor- 
tices in-situ density distributions 0], and the con- 
densate fraction [J, 0] established a critical polariza- 
tion P c = (iVj - N^/Ntot = 0.70(3) for the superfluid 
phase to exist. These results were obtained in a har- 
monic trap with cylindrical symmetry (uj x = uj y — oluj; 
u> z = lo), with aspect ratio a ~ 5, and total particle 
numbers N tot = +N i ~ 10 6 - 10 7 . The Rice experi- 
ment 0, [f| also observed phase separation, with a fully- 
paired core surrounded by normal majority fermions, 
but with a sharp phase boundary and extremely thin 
partially-polarized shell for low temperatures, and the 
unpolarized core exists to high polarization P c > 0.9. 
These experiments are in a highly elongated trap, with 
aspect ratio a ~ 35 — 45, and for lower A^ot ~ 10 5 . For 
these conditions, the core deformation and the double- 
peak structure in the axial density imply a breakdown of 
the local-density approximation (LDA) in Refs. 0, HJ . 

The critical polarization is influenced by the energy of 
the competing normal polarized phase. For large asym- 
metries, this is governed by the energy of a spin-down 
fermion interacting resonantly with a spin-up Fermi sea, 
which is a universal function of a and N . We study 
the energy of this spin-down fermion, the so-called Fermi 
polaron, in anisotropic traps for different particle num- 
bers and show that the experimental differences can be 
understood partially based on our microscopic results. 
This strongly-interacting Fermi polaron provides insights 
to problems in condensed-matter systems, with lower di- 



mensions playing the role of trap-geometry effects, as well 
as to nuclear physics, where neutron-rich nuclei exhibit 
neutron skins 0], with neutron/proton densities similar 
to the spin densities in resonantly-interacting cold atoms. 

Uniform system. The polaron energy E was calcu- 
lated variationally for the uniform system including one- 
particle-one-hole excitations (lplh) and estimated 
in Ref. [HI]. This leads to a Schwinger-Dyson equation, 
G^ 1 (E) = 0, or diagrammatically: 
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where G (Go) is the full (noninteracting) spin-down 
propagator with momentum p = 0. For large S-wave 
scattering lengths, l/a s — 0, the energy is universal, 
E = hi = r/p F /(2m), with Fermi momentum pp. The 
self-consistent solution to Eq. ([I]) yields r\ = —0.607 [9j. 

This energy gain constrains the equation of state for 
large asymmetries, and thus the existence of partially- 
polarized phases and the critical polarization: The vari- 
ational 77 is lower than the maximal stress (from /ij — 
fit ^ 2 A) for stability of the superfluid phase, Ml/MT ^ 
—0.09(3), and this requires the existence of at least one 
nontrivial partially-polarized phase in the uniform sys- 
tem d, [n|. In LDA with Eq. |{I3|), r\ = -0.607 leads 
to P c = 0.74 and a critical density ratio x c = ni/n-\ = 
0.47 (HI, which are in good agreement with P c of the 
MIT experiment [H, [H, 3 and with a tomography mea- 
surement of x c w 0.47 [TJ. Finally, the variational lplh 
77 value agrees very well with Monte-Carlo (MC) re- 
sults [13, 13 . 3, EH , and 2p2h contributions were shown 
to be small [161 ]. 

Basic formalism. The strongly-interacting Fermi gas 
in a harmonic-oscillator trap is given by the Hamiltonian 

11,(7 
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where e n = otu)(n x + n y + 1) + uj{n z + 1/2) are harmonic 
oscillator energies (K = 1). The operator a n ,cr annihilates 
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a particle with spin a =| , J, in a state with quantum num- 
bers n = (n x ,n y , n z ). We use a contact interaction regu- 
lated by separable cutoff functions in momentum space, 

(p| W> = C(A) e-^'^ with C(A) = 3^3- , 

(3) 

where p, p' are incoming/outgoing relative momenta, m 
is the fermion mass and A a momentum cutoff. In this 
case, the harmonic-oscillator matrix elements can be ex- 
pressed as a sum over separable functions F(ja.\, ri2, S), 

(m, n 2 |V> 3 , n 4 ) = C(A) ^ F(m, n 2 , S)F(n 3 , n 4 , S) , 



with center-of-mass quantum numbers S, -F(n l7 n 2 , S) 
dimcnsionlcss function F is given by 



(4) 
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f(ni,Si,n 2i ), (5) 



where the relative quantum numbers = ni i + n 2i — 
have to be even and positive, and one has for n 2i ^ Tij 



f(n l ,S i ,n 2l ) 



rii 
n 2i 



2 Fi (-7l 2j , — -Si, 1 - 7l 2 , +Tii, -1) , 

(6) 



with hypergeometric function 2 i*i, and for n 2i > rij 



x 2 Fi(-rii, n 2i - Hi - 



Si 

n 2z - n-i 
Si, l + n 2i - ri,, -1) 



(7) 



Polaron energy. Following the variational Ansatz of 
Refs. we calculate the energy _B of the spin-down 

fermion, including lplh excitations in the wave function, 

M=6)|fi>+ ^ m , h , p |m,h,p>, (8) 

m.h.p 

where denotes the Fermi sea with the spin-down par- 
ticle in the n = level 20], and |m, h, p) consist of a 
spin-up fermion in h excited to a level p above the Fermi 
energy Ep, and the spin-down particle occupies the level 
m. Therefore, the sum over h is restricted to occupied 
states, whereas p is over unoccupied states above e-p. 

Minimizing (ip\H\ip) with respect to </>o, </3 m ,h,p, we find 
the self-consistent equation for E in anisotropic traps, 

E-e = 

]T ]TF(0,h,S) [Af- 1 (e Pj £; + eh )] SiIi iJ'(0,h,L), 

(9) 




FIG. 1: (Color online) Upper panel: Local Fermi energy 
Ef(1,N) at the center of an isotropic trap divided by the 
large- TV expression £f,oo(1, TV) as a function of spin-up parti- 
cle number TV. Lower panel: Energy n(l,N) for an isotropic 
trap, with fits to our numerical results (see text). The hori- 
zontal line represents r\ = —0.607 for the uniform system 0. 



where E is measured from the energy of the Fermi sea, 
in weak coupling E rs eo, and the matrix M is given by 



M(e F ,E + e h 



S,L 



1 



C(A) 



- D(a,AE) 



ds.L 



^ ^ E + e h -(e p +e m ) 



(10) 



Here AE = a(S x + S y + 2) + S z + 1 - (E + e h )/w and 



D(a, AE) is identical to the last term of Eq. (fl~0"|) with 
unrestricted sum over p and S = L. For an isotropic 
trap, D(l, AE) has the simple analytical form 



D(1,AE) 



lA 



2(2^) 3 / 2 



raw 
"2V 



3/2 



V^r(A£/2) 



loT((AE -l)/2) 

(11) 

The cancellation of the first term in Eq. (fTl"]) with the 
cutoff in the 1/C (A) term in Eq. (jTUJ) demonstrates that 
E is cutoff independent for large A. We have verified that 
this is the case for all studied a and use A > 10 4 y^mw/2. 
Moreover, we have found numerically that the diagonal 
matrix elements of M(ef,E + eh) depend only on the 
ccntcr-of-mass excitation a(S x + S y ) + S z . 

For large scattering lengths, l/a s = 0, the energy is 
a universal function of the aspect ratio and the spin-up 
particle number TV — TVj, and we generalize the scaling 
for the uniform system to anisotropic traps, 



E = r)(a, TV) E F (a, TV) . 



(12) 
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FIG. 2: (Color online) Energy n{a, TV) as a function of TV" 
for various aspect ratios a, compared to 77 = —0.607 for the 
uniform system [9( (horizontal lines). The fits are discussed 
in the text and shown separately in the lower panel. 



where E F (a, TV) = (67r 2 n T (0)) 2/3 /(2m) is the local Fermi 
energy of spin-up particles at the center of the trap. 
In the upper panel of Fig. [TJ we show E F (1,N) for 
an isotropic trap divided by the large-TV expression 
E F<00 (1,N) ^uj{6N)^ 3 . The points are for alternating 
odd-even values of the Fermi level n F , which defines the 
Fermi energy e F = u>(an F + (2a + l)/2). The local Fermi 
energy approaches E Fi00 (l,N) from above (below) for 
odd (even) n F . With increasing a, this effect decreases 
and the envelopes approach the large- TV result faster. 

Results. Using Eq. (fT2"|) . we solve Eq. © iteratively 
for r](a,N), with a numerical precision better than 1%. 
To this end, we take the matrix M(e F ,E + eh) to be 
diagonal. This is correct in the large-TV limit, and we 
have checked numerically that the off-diagonal matrix el- 
ements are considerably smaller than the diagonal ones 
for all studied values of TV. In the lower panel of Fig.[IJ we 
show 77(1, TV) for an isotropic trap as a function of TV. The 
effects due to finite particle numbers and confinement of 
the trap are clearly present: The odd-even systematics 
seen in the local Fermi energy E F (1,N) is small com- 
pared to the decrease of 77(1, TV) with particle number. 
Therefore, the decrease is not due to the change in the 
local Fermi energy. For one s pin- up fermion, the exact 
ground-state energy in a trap [17j is E(l, 1) = ui/2 > 0, 
thus 77(1,1) = (3671-)- 1 / 3 = 0.207 [2l|. With increasing 
TV, 77(1, TV) decreases and saturates. Using the Ansatz, 
77(a,TV) = a(a)(l + 6(a)TV- c ( Q )), wc fit our numerical re- 
sults for odd (even) n F separately and find a(l) —0.61 
and c(l) « 0.34 (0.32). This is in very good agreement 



FIG. 3: (Color online) Upper panel: Critical polarization 
P c (a,TVt t) as a function of aspect ratio a for TVtot = 10 4 
(upper) and TV to t = 10 (lower set of curves). Lower panel: 
Critical density ratio x c (a, N to t) for TVtot = 10 4 (lower) and 
TVtot = 10 5 (upper set of curves). Results are shown for two 
approximations to the quasiparticle spectrum and interaction. 



with 77 = —0.607 for the uniform system [9( and natural 
large-TV corrections of 1/E F<00 (1, TV) ~ TV" 1 / 3 . There- 
fore, ~ 10% changes of 77 are natural for TV ~ 10 4 . 

In Fig. [51 we show the dependence of rj(a, TV) on trap 
geometry, for various aspect ratios from a = 1 to a = 
35, as a function of the spin-up particle number. The 
single-particle energy depends significantly on the aspect 
ratio, while the odd-even Fermi level effect decreases with 
increasing TV, a and is negligible for a > 10. For fixed 
TV, 77(0;, TV) increases with increasing a. In addition, for 
larger aspect ratios, the dependence on TV is stronger. 
For each a, we fit our combined results (including odd 
and even n F ) with the power-law Ansatz and show the 
fits in Fig. [2] We find a(a) ss —0.61(1), consistent with 
the uniform result for all studied aspect ratios, and c(a) 
ranges from c(l) ~ 0.36 to c(35) « 0.31. 

Critical polarization. We now explore the impact of 
the calculated finite-size and confinement (trap) effects 
on the phase structure. We consider an unpolarized su- 
perfluid phase and a partially-polarized normal Fermi liq- 
uid. Following Ref. [l9(, the free energy is given by 

E tot = 2 [e s K(r)) + V(r) - f i s ] n s (r) dr 

J\r\<Ra 



ife<|r|<i?. T 



£ n (ir(r))n T (r) + U(r)(n i (r) +n T (r)) 



/J T n T (r) 77|(r)] dr., 



(13) 



where x — njn^ ^ 1, Rg = a(R 2 x +Ry) + R 2 z defines the 
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boundary of the superfluid phase, and the excess spin-up 
density vanishes at i?| . As discussed, the LDA of Eq. (fl3|) 
breaks down for the Rice experiment 0, 0] • We only use 
this here to explore the impact of rj(a, N) on the critical 
polarization. In a full density-functional calculation, this 
can also be combined with surface tension or gradient 
terms. For the uniform system at unitarity, the energy 
density of the superfluid e s and of the partially-polarized 
normal Fermi liquid e n are given by [191 ] 



Ms 



^5 



3 (67r 2 n s ) 2 /3 



2m 



and e n (x) = — 
5 



3(67r 2 n T ) 2 / 3 



2 to 



e{x), 
(14) 

with superfluid density n s and universal energy £ = 0.42 
of the symmetric system 13|, |15[ , which is consistent with 
£ = 0.46±0.05ofRef. @. Assuming i<1, the energy of 
adding spin-down fermions to the normal phase is deter- 
mined by rj(a,N), with corrections due to a spin-down 
quasiparticle effective mass to* and due to quasiparticle 
interactions B 1191: 



e(x) 



1 + — rj(a, N) x H a 

3 m 



5/3 



Bx' 



(15) 



We take r](a, N) from Fig. [21 but for simplicity consider 
two cases for the quasiparticle spectrum: m*/m = 1, 
B = 0, as well as the MC values to* /to — 1.09, B = 
0.14 [15j, which show these are corrections to the leading 
effects from r\. This however does not include the effects 
of Fermi statistics of the minority particles on r\. 

The critical polarization P c (a, AT tot ) is obtained, when 
the phase boundary reaches the trap center R$ — > 0. In 
chemical equilibrium, fig = + Mt)/2j the ground state 
of the system is determined by requiring that the energy 
functional, Eq. (|13p . is stationary with respect to vari- 
ations of the densities and of the phase boundary i?s, 
so that the pressure between the two phases is equal: 
2n1{de s / dn s ) — n 2 (9e n /3ri|) + n^n^de^/ dn{). This 
leads to an equation for the critical density ratio at the 
center of the trap, x c (a, N to t) = x(Rs = 0,a,N to t) [H3]: 
e(x c ) + |(1- x e ) d x e{x c ) - (2£) 3/5 e(x c ) 2 / 5 = 0. Given 
iVtot and x c , the spin-up/spin-down densities and particle 
numbers are determined from the variation of E to t- 

In Fig. [3l we show the dependence of the critical 
polarization P c (a, N to t) and the critical density ratio 
x c (a, N to t) as a function of aspect ratio, for total particle 
numbers Ntot = 10 4 and AT to t = 10 5 , where the experi- 
mental differences from the uniform system exists. For 
a = 1, Ntot — 10 7 and the MC m* , B values, we reach 
the uniform system P c — 0.74 and x c = 0.47. For fixed 
Ntot = 10 4 and increasing a from 1 to 35, P c increases 
from 0.82 to 0.89 and x c decreases from 0.40 to 0.31 (for 
the MC to*, B values). For fixed a — 35, P c increases 
with decreasing N tot = 10 5 , 10 4 , 10 3 from 0.82, 0.89, 0.96 
and x c decreases from 0.39, 0.31, 0.18 (for the MC to*, 
B values; results for iVtot = 10 3 not shown in Fig. [3j). In 
addition, we show in Fig. [3] the dependence on the quasi- 
particle spectrum (through m*) and interaction B. For 



given a and N to t , Pc is larger and x c smaller for the MC 
to*, B values, compared to to*/to = 1, B = 0, but as ex- 
pected, the uncertainty due to to*, B is smaller than the 
variation of P c and x c with a, N to t- This dependence also 
becomes weaker with increasing a and decreasing iV tot . 

In summary, for lower particle numbers and more elon- 
gated traps, the energy of the normal polarized phase in- 
creases and the superfluid extends to larger population 
imbalances. This provides a microscopic understand- 
ing of the MIT-Rice differences due to the dependence 
of the polaron energy on the particle number and the 
trap geometry. Finite-size effects are stronger in highly- 
elongated systems as the dimensionality of the problem is 
continuously reduced with increasong aspect ratio. The 
N + 1-body problem is a natural first step towards gen- 
eral asymmetries and towards contributions to the total 
energy beyond T)(a,N). In addition, effects from a full 
density-functional calculation need to be studied. 
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The restriction to n = (justified for large N) enables 
our still involved numerical solution. 
The exact two-body ground-state energy [TtJ is repro- 
duced, if we generalize the first term in Eq. ([8} to include 
a sum over the spin-down particle in level n, <p n \n). 



